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COMPUTER EXPERIMENTS ON THE STRUCTURE 
AND DYNAMICS OF SPIRAL GALAXIES 


By Frank Hohl 
Langley Research Center 

SUMMARY 

The evolution of an initially balanced rotating disk of stars with an initial velocity 
dispersion given by Toomre's local criterion is investigated by means of a computer 
model for isolated disks of stars. It is found that the disk is unstable against very large- 
scale modes. After about two rotations the central portion of the disk tends to assume 
a bar-shaped structure. A stable axisymmetric disk with a velocity dispersion much 
larger than that given by Toomre's criterion is generated. The final mass distribution 
for the disk gives a high-density central core and a disk population of stars that is closely 
approximated by an exponential variation. 

Various methods and rates of cooling the hot axisymmetric disks were investigated. 
It was found that the cooling resulted in the development of two-arm spiral structures 
which persisted as long as the cooling continued. 

An experiment was performed to induce spiral structure in a galaxy by means of the 
close passage of a companion galaxy. Parameters similar to those expected for M51 and 
its companion were used. It was found that because of the rather high velocity dispersion 
of the disturbed disk galaxy, only a weak two-arm spiral structure appeared as the result 
of the passage. 

The evolution of a uniformly rotating disk galaxy which is a stationary solution of 
the collisionless Boltzmann equation is investigated for various values of the initial rms 
velocity dispersion. It is found that the disk becomes stable at a value of the velocity 
dispersion predicted by theory. 


INTRODUCTION 

With the introduction of computer models for thin self- gravitating stellar systems, 
the field of experimental stellar dynamics is beginning to provide fresh insights into the 
structures of spiral galaxies. Two computer models for self-consistent disk galaxies 



have recently been described, and to a certain extent demonstrated, by Miller and 
Prendergast (ref. 1) and by Hohl and Hockney (ref. 2.). 

To investigate the development of spiral structure, the model was modified by 
Hohl (refs. 3 and 4) to include a fixed central force similar to the Schmidt model of the 
Galaxy. It was then found that spiral structure persisted for more than eight rotations. 
Miller, Prendergast, and Quirk (ref. 5) have also modified their model to include a dissi- 
pative component resembling gas. They found that some very interesting spiral patterns 
developed, especially in the "gas" component, and that such patterns remained for about 
three galactic rotations. Such a combined system of stars and gas seems more realistic 
for the study of the spiral phenomena than a model composed only of stars. 

However, a basic understanding of the structure of galaxies clearly requires a 
knowledge of how purely stellar disks would behave. Accordingly, the present report 
will be concerned primarily with the dynamical evolution of isolated disks consisting of 
stars alone. The computer model used for the calculations is that developed by Hohl and 
Hockney (ref. 2). From previous work (refs. 3 and 6 to 9) it is, of course, already clear 
that disks of stars with velocity dispersions less than those estimated by Toomre (ref. 6) 
to be locally stabilizing are Violently unstable. For this reason the present paper reports 
primarily on disks with initial velocity dispersions at least equal to those given by 
Toomre's local criterion. In a previous paper (ref. 9) the authors were overenthusias- 
tic in calling such disks stable . It would have been more appropriate to say that only the 
fast-growing small-scale gravitational instabilities are avoided when Toomre's criterion 
is satisfied. As will be seen herein, such disks remain susceptible to slower large-scale 
instabilities, primarily of the bar -making type. 

Thus, disks of stars are considerably more difficult to stabilize than indicated by 
local analyses. It is, of course, only with the advent of computer experiments such as 
those described here that it has become possible to determine the behavior of disks of 
stars in the large. Before such experiments, there existed no really adequate theory 
describing the overall dynamics of disks of stars. 

Freeman (ref. 10) has summarized the observational evidence that most spiral and 
SO galaxies have a dense spheroidal central component and a disk component of stars 
described by an exponential light distribution. One of the questions the present report 
tries to answer is whether a disk of stars will evolve to a similar final state. 

In published work (refs. 1, 2, 3, 5, and 9) on the simulation of disk galaxies, the ini- 
tial conditions chosen were not solutions of the time-independent collisionless Boltzmann 
equation. Even though some of the initial conditions (refs. 3 and 9) were for initially bal- 
anced disks, these initial conditions did not correspond to stationary disks. One excep- 
tion to this is the "cold" balanced disk, which, however, is violently unstable. In the 
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present paper the evolution of initially stationary disks for various temperatures is 
investigated. 

The computer model used for the calculations is described in the appendix. 

How striking the resemblance of computer- simulated disk galaxies is to actual 
galaxies is shown in figure 1. In that figure four computer- generated galaxies are com- 
pared with photographs of actual galaxies. 

SYMBOLS 


E total energy 

f distribution function 

G gravitational constant 

H Green's function 

J angular momentum 

K gravitational field, V<£ 

M 0 mass of Sun 

m mass 

N dimension of array used in potential calculations 

n = N/2 
Q = tfr/o^min 

R radius of disk galaxy 

r,0,z cylindrical coordinates 

T q rotational period of cold balanced disk 

t time 
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V 


rotational velocity 


v 

x,y,z 

6(z) 

K 

V 

a 

a r,min 

<P 

co 

o 

Subscripts: 

i,j 

k ,1 
max 
r,e 
x,y 

Notation: 


velocity 

Cartesian coordinates 
Dirac delta function 
epicyclic frequency 
surface mass density 
velocity dispersion 

velocity dispersion defined by equation (3) 
gravitational potential 
angular velocity 

angular velocity of cold balanced disk 

summation indices 
summation indices 
maximum 

radial and azimuthal component 
summation indices 

Fourier transformed quantity 
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EVOLUTION OF AN INITIALLY BALANCED UNIFORMLY ROTATING DISK 


In the following sections the evolution of a balanced uniformly rotating disk of stars 
with an initial velocity dispersion given by Toomre's local criterion (ref. 6) is investi- 
gated. Disks with other initial mass distributions have previously been investigated in 
less detail (ref. 3). 


Initial Conditions 

The present report describes in detail the evolution of a balanced uniformly rotating 
disk of stars with an initial velocity dispersion given by Toomre's local criterion. The 
surface mass density of the uniformly rotating disk is given by 


ju(r) = jn(O). 


1 


r2_ 

R 2 


( 1 ) 


where r is the radial coordinate, R is the radius of the disk, and /i(0) is the central 
surface mass density. The uniform angular velocity required to balance the cold (zero- 
velocity-dispersion) disk is 



Gfi(O) 

2R 


( 2 ) 


where G is the gravitational constant. 

Toomre’s local criterion (ref. 6) for suppression of all axisymmetric instabilities 
requires a minimum radial velocity dispersion (Gaussian velocity distribution) given by 


a 

r,min 


3.36 


Gjii. 


(3) 


where a r 2 is the second moment of the distribution function, p is the local value of 
the density, and k is the local value of the epicyclic frequency. The epicyclic frequency 
is defined by 


/c 2 = 




(4) 


where K r is the radial component of the gravitational field. For the uniformly rotating 
disk K r = a> 0 2 r; thus, k = 2 oj q and 
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( 5 ) 


cr ■ = 0.341Ra>„ 

r,min o 



To determine whether Toomre's criterion will stabilize the disk, the initial (Gaussian) 
velocity dispersion is taken as 


°r( r ) °r,min 

for the radial component and 


( 6 ) 


a 0 (r) = 


fl(r) 

2w Q (r) 


°r,min 


(7) 


for the azimuthal component. For the uniformly rotating disk, equations (6) and (7) 
become 


o- r (r) = cr fl (r) = 0.341Ra) o ^/l - (8) 

Because of the added velocity dispersion the initial angular velocity of the stars u> for 
a balanced disk is lower than w 0 . By summing the forces acting on a small surface ele- 
ment of the disk, the following equation is obtained: 

w2 = W ° 2 + rfer ^ (r)cT r 2 ( r >] + r-2[a r 2(r) - <r 0 2 (r)] (9) 

For the uniformly rotating disk, equation (9) gives w = 0.809u> o , which results in a cir- 
cular velocity given by 

V Q = 0.809rw o (10) 

The initial conditions do not exactly satisfy the collisionless Boltzmann equation for a 
stationary stellar disk. An initial condition with a Gaussian velocity distribution was 
preferred to the stationary initial state corresponding to the density given by equation (1), 
because Toomre's criterion is strictly valid only for Gaussian velocity distributions. The 
stationary state for the uniformly rotating disk (ref. 11) has a peculiar non-Gaussian veloc- 
ity distribution such that the density in phase space is not a decreasing function of epicy- 
clic amplitude; hence overstabilities may be present (ref. 12). With the exception of the 
uniformly rotating disk, there is a rather unfortunate lack of exact, hot equilibrium models 
for disks of stars. 
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The initial distribution of stars is generated by means of a pseudorandom -number 
generator. The initial circular velocity is given by equation (10), and superposed on the 
circular velocity is a velocity dispersion with a Gaussian distribution corresponding to 
°r,min as S iven by equation (8). 


Axisymmetric Evolution 

Since the initial state chosen for the model is not an equilibrium state, the axisym- 
metric behavior of the disk is first investigated to determine whether any axisymmetric 
instabilities or other misbehavior might be present. This is done by considering the evo- 
lution of the initial disk under a gravitational field that is constrained to remain purely 
radial. This gravitational field is obtained by averaging at each time step the fields com- 
puted in the ±x-directions and the ±y-directions. The results are shown in figure 2 for 
eight rotations. Initially the disk has an arbitrarily chosen radius of 8 kpc. The square 
border enclosing the disk is then at x = ±19 kpc and y = ±19 kpc and covers only part 
of the 128 x 128 square array of cells used for the gravitational-field calculations. Each 
of the 100 000 stars in the disk can be thought of as having a mass of 0.84 x 10®M o . The 
time shown in figure 2 and in subsequent figures is in units of the rotational period of the 
cold balanced disk; each such unit equals 1.5 x 10** years. 

The disk shown in figure 2 remains, of course, axisymmetric. After a few small 
pulsations it settles down to an essentially steady state. When the total kinetic energy of 
the disk is plotted as a function of time, the initial amplitude of the oscillations in the 
total kinetic energy is about 8 percent of the mean value. These oscillations in the total 
kinetic energy are quickly damped; after four rotations the amplitude of the oscillations 
is less than 1 percent of the mean value. Figure 3 shows the evolution of the density. 
Again it can be seen that the system quickly reaches a steady state . During the evolution 
of the disk, the value of ° v / a v m i n remains near unity for r less than 5 kpc. These 
results show that any troubles the disk may have are basically nonaxisymmetric in nature. 

The disk in figure 2 at t = 8 should now be very close to a stationary state of the 
collisionless Boltzmann equation. If the evolution of the disk is continued without the con- 
straint of keeping the gravitational field radial, the evolution shown in figure 4 is obtained. 
In less than two additional rotations the disk assumes a bar- shaped structure - an indi- 
cation that the newly made equilibrium is unstable in a nonaxisymmetric sense. 

Fully Nonaxisymmetric Evolution 

These instabilities are now studied by returning to the smooth but slightly nonequi- 
librium disk that was the starting point for figure 2. Those initial conditions, by the way, 
are such that the virial theorm is satisfied — that is, the negative of the calculated poten- 
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tial energy equals twice the kinetic energy with an accuracy of about 0.5 percent. This 
small discrepancy is meaningless because of the approximate method used for calculating 
the potential energy (ref. 2). Also, the maximum change in the calculated total energy of 
the system during its evolution is near 1 percent, and the variation in the total angular 
momentum is even less. 

The totally unconstrained evolution of this disk is displayed in figure 5. It shows 
that the disk is indeed stabilized against small-scale disturbances which would completely 
disrupt a cold disk in less than one rotation (ref. 3). However, it is again found that the 
disk is not stabilized against relatively slow-growing large-scale disturbances which 
cause the system to assume a very pronounced bar-shaped structure after two rotations. 
After about four rotations the disk population of stars constitutes a nearly axisymmetric 
distribution surrounding a dense central oval or bar-shaped core. The oval contains 
nearly two-thirds of the total mass of the disk and has an axis ratio of about 3:2. After 
about four rotations there is little change in the structure of the disk, and the bar rotates 
with a constant period of 2.25 times the period of the cold balanced disk. 

To obtain more quantitative information on the structure and evolution of the bar- 
shaped central part of the system, a number of annular rings 1 kpc in width and concen- 
tric with the disk were superposed on the disk shown in figure 5. Each annular ring was 
divided into 100 equal segments, and the number of stars in each segment was determined. 
The results are shown in figure 6 for annuli centered at the three radii shown. The fluc- 
tuations caused by generating the initial positions of the stars by means of a random- 
number generator can be seen for the curves at t = 0. At this time there are no stars 
in the annulus at r = 10 kpc because the initial disk has a radius of only 8 kpc. After 
one rotation the bar structure begins to form and stars are ejected beyond the 10 kpc 
radius. At t = 2.0 the bar is very pronounced, and after this time it begins to broaden. 
After eight rotations the disk has assumed a steady state, as indicated by the result that 
neither the amplitude nor the structure of the rotating bar changes any longer. 

To picture the radial variation of parameters describing the disk, the disk was 
divided into a number of concentric rings each 0.5 kpc in width. The radial dependence 
of various parameters averaged azimuthally over each ring was then obtained. Figure 7 
shows the radial variation of the rms radial-velocity dispersion obtained in that manner. 

At t = 0 the radial velocity dispersion is described by equation (5). As high-velocity 
stars are ejected radially outward, the velocity dispersion becomes large for large radii. 
After about eight rotations the velocity dispersion remains essentially unchanged. In 
interpreting the results shown in figure 7 and in subsequent figures, the nonaxisymmetric 
shape of the disk for small radii (r < 8 kpc) should be kept in mind. A better indication 
of how hot the stars in the disk are can be obtained from a plot of the azimuthally aver- 
aged value of 
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Q 


g r (r) 

°r,min 

(The value of cr r min = 3.36 G/i//c is obtained by using n as given in fig. 9 and by 
obtaining k from the slope of K r .j The radial dependence of Q is shown in figure 8. 
Initially Q is given a value near 1. However, as the disk evolves, Q increases rapidly 
in the outer regions of the disk. Figure 8 is typical of the high velocity dispersion in the 
outer parts of the disk for computer- generated galaxies. For a number of other initial 
conditions (ref. 3), typical values of Q were found to be between 2 and 4. 

The evolution of the surface mass density is shown in figure 9. It can be seen from 
this figure that the central mass density increases by about a factor of 4 during the 
11.6 rotations shown. After eight rotations the radial variation in density changes only 
very little. To determine how closely the final density variation corresponds to an expo- 
nential variation, the density for t = 8 and t = 11.6 rotations is plotted in figure 10 on 
a semilog scale. This distribution of the disk population of stars (outside r = 8 kpc) is 
closely approximated by an exponential variation with a scale length of 8.6 kpc. In the 
high-density central oval the density also closely follows an exponential law, but with the 
much shorter scale length of 1.44 kpc. Figure 10 also illustrates how little the structure 
of the disk changes during the last 3.6 rotations. The density functions shown in figures 9 
and 10 represent the azimuthally averaged value. In figure 11 is presented the unaveraged 
density at t = 8 along the bar and across the bar. As expected, the density falls off 
faster across the bar than it does along the bar. For radii greater than about 10 kpc the 
densities in the two directions are essentially the same. 

GENERATION OF A STABLE AXISYMMETRIC DISK 

It was found that disks of stars with initial conditions generated according to analyt- 
ical expressions such as equations (1) to (10) are generally unstable and finally assume a 
steady state with a central bar-shaped structure. The evolution of such warm disks is 
found to be very similar, irrespective of whether the initial mass distribution is Gaussian, 
exponential, or some other distribution (ref. 3). Also, the stationary state finally reached 
by the disks results in a very hot population of stars in the outer part of the disk, as is 
shown in figure 8. In order to generate an axisymmetric stable disk, the final distribution 
of stars (at t = 11.6) of the disk shown in figure 5 was used as an initial condition after 
symmetrizing out the bar structure. All stars kept the same (radial and tangential) veloc- 
ity components and radius, but they were redistributed randomly in longitude. This new 
disk has the same radial variation of its different parameters as given (at t = 11.6) in 
figures 5, 7, 8, 9, and 10. 
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The reason for using this highly artificial procedure of "symmetrizing” was mainly 
to show that the bar structure, although here persistent, is not necessarily very strong or 
"deep." Figure 12 depicts the evolution of the new disk for six rotations. The disk is 
now axisymmetric and stable; all parameters remain essentially constant during the evo- 
lution. For example, the variation of the density given in figure 13 shows no change 
between two and six rotations. Note also that the "experimentally" obtained density vari- 
ation is closely approximated by the sum of the exponentials 

M(r) = 4.2 x 109 exp(- ^ + 3.6 x 10? exp(- ^ (11) 

Since the disk is axisymmetric and stable, the variation of its rms velocities near a fixed 
radius should provide information on collisional or thermalization effects that may be pres- 
ent in the model. In figure 14 the rms velocities of stars in a 0.5-kpc-wide ring centered 
at the four radii shown are plotted as a function of time. It can be seen that there is no 
systematic increase (or decrease) in the rms velocities for the six rotations shown. 

Figure 15 shows the radial dependence of the mean circular velocity of the stars 
(Vg), of oj ~ ^K r /r , and of rw. It is interesting to note that the azimuthally averaged 
radial dependence of the parameters of the disk did not change after the central bar struc- 
ture was eliminated by symmetrizing. The values of Q for the axisymmetric disk are 
the same as shown in figure 8 at t = 11.6 and are rather high. For some of the disks 
investigated, final values of Q near 2 were obtained. However, for most of the disks 
investigated (ref. 3) the final value of Q is near 4. 

COOLING OF THE AXISYMMETRIC DISK 

The stable axisymmetric disk described in the previous section is rather hot, espe- 
cially in the outer regions of the disk. It is therefore of interest to determine the effects 
of cooling the disk. Various methods of cooling were tried. For example, during each 
rotation a certain percentage (from 5 to 30 percent) of the stars, chosen at random, had 
a portion of their noncircular velocities removed in proportion to their radii. Thus stars 
near the center kept nearly all their random motion whereas stars far away from the 
center were placed in more nearly circular orbits. Another method, somewhat more 
analogous to the "gas collisions" of Miller et al. (ref. 5) was simply to place a certain 
percentage of the stars during each rotation into purely circular orbits. All methods of 
cooling which were tried gave essentially the same results. 

Figure 16 shows the effects of cooling when for each rotation 10 percent of the stars 
are selected at random and are placed in circular orbits. The initial disk is that shown 
in figure 12 at t = 6. During the first three rotations the appearance of the disk changes 
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little. However, at t = 3.1 there is already a pronounced bar structure at r = 2.5 kpc. 
After t = 3.1 a two-arm structure appears, which during the following 1.5 rotations dis- 
plays a theta- shaped structure. After t = 4.5 the spiral structure opens up and rotates 
with a fairly constant speed of about 10 km s~l kpc - ^. This compares with a mean circu- 
lar velocity of the stars of about 15 km s - l kpc~l at r = 10 kpc. These results tend to 
support the density-wave theory of Lin (ref. 13). At t = 6 the value of Q for the disk 
was between 2 and 3. Any further cooling only caused the collective instabilities to heat 
up the disk as fast as it was being cooled. This effect can be seen in figure 17 where the 
radial velocities of the stars are plotted as a function of radius. At t = 2 the effect of 
placing a number of stars in purely circular orbits is clearly visible by the concentration 
of stars with zero radial velocity. At t = 4 an increased velocity dispersion corre- 
sponding to the condensation of stars in the theta-shaped spiral arms is clearly visible. 
Also, the effects of heating by means of collective instabilities prevent a reduction of the 
radial velocity dispersion after t = 6. Figure 18 displays the azimuthal variation of the 
number density at three radii. These results show that in addition to the density-wave- 
like spiral structure of the outer disk, the galaxy also develops a strong barlike structure 
of the central core. 

The results of removing a portion of the stars’ random velocity in proportion to the 
star's radial distance is shown in figure 19. Thirty percent of the stars are affected 
during each rotation. The results are slightly different from those of figure 16, primar- 
ily because the cooling rate is now larger. In fact cooling proceeds now so rapidly that 
at t = 11 clusters of stars have condensed and are circling the central core as satellites. 

EFFECT OF SATELLITE GALAXY 

Alar Toomre of the Massachusetts Institute of Technology has studied a number of 
spiral galaxies and found that many have nearby companion galaxies. To determine the 
effect of a close passage, a stable axisymmetric disk galaxy similar to that shown in fig- 
ure 12 was perturbed by the passage of a companion galaxy having one-fourth the mass of 
the primary galaxy. The parameters for M51 were chosen for the calculations. The orbit 
was direct, the inclination of the companion's orbit plane to the plane of the primary gal- 
axy was 75°, the angle between the node and periapsis was 20°, and periapsis was 25 kpc. 
The effect of the passage on the structure of the primary galaxy is shown in figure 20. 

The time is given in rotational periods of the primary galaxy; t = 0 corresponds to the 
time at periapsis. These results indicate that only a weak two- arm spiral structure 
developed. This was to be expected since the velocity dispersion of the primary galaxy 
is rather high, as shown in figure 21. The value of Q ranges from about 2 in the cen- 
tral region to 6 in the outer disk. A more detailed view of the velocity distribution is 
given in figure 22, where the radial velocity distribution for the stars in various annular 


11 



rings is given. The effect of the large random velocities on the formation of spiral struc- 
ture can be determined by plotting the distribution of stars in various velocity intervals. 
This is done in figure 23. As can be seen, the spiral structure is quite pronounced for the 
lower velocities, but can hardly be detected for the higher velocities. This was to be 
expected since the large random velocities have a dispersive effect on the formation of 
spiral structure. It thus appears that a much cooler primary galaxy is needed to study 
the effects of a close passage. 

SELF-CONSISTENT STATIONARY DISKS 


Presently, there are no really "good" stationary solutions of the collisionless 
Boltzmann equation available for disk galaxies. For certain mass distributions analyti- 
cal stationary solutions can be obtained; however, these solutions have certain undesirable 
features (e.g., the distribution function becomes singular at the edge of the disk) which 
make them unlike any "real" galaxy. Nevertheless, there are analytical predictions as 
to the stability and behavior of such disks. Thus, by investigating such stationary states, 
it is possible to determine whether the analytical methods used to study the behavior of 
disk galaxies yield the correct predictions. 


Theory 

A time-independent, axisymmetric self-consistent disk of stars is described by the 
collisionless Boltzmann equation 


, 8f. V 

r 9r r 8v r 


v r v fl at 
r 8v 0 



( 12 ) 


where v r and v^ are the radial and azimuthal velocity components and f/r,v r ,vA is 
the distribution function defined so that 


dm = f^r,v r ,v 0 ^dv r dv 0 r dr do (13) 

corresponds to the mass in a phase-space element within the surface element r dr dd. 
The gravitational field K r equals where the potential <p( r,z) is obtained from the 
Poisson equation 

v2 (f) = 47rG p(r) 6(z) (14) 

where G is the gravitational constant and the density ji (r) is given by 
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» I < r) =Jj''( r > v r> v e) dv r dv 0 


(15) 


According to Jeans' theorem (refs. 14 and 15), any function of the form 

f(r,v r ,v 0 ) = f(E,J) (16) 

is a solution of equation (12), where the energy E and the angular momentum J are 
given by 

E = |( v r 2 + V) + * 

and 

J = rv 0 

An axisymmetric mass-density variation for which the corresponding distribution 
function can be obtained is that for the uniformly rotating disk 

where p. (0) is the central mass density and R is the radius of the disk. The gravita- 
tional potential inside the disk is (ref. 3) 

0(r) = j w o 2 r2 - u> 0 2 R 2 (17) 


where co 0 = 77|/Gp(0)/2R is the uniform angular velocity required to balance the cold 
(zero-velocity-dispersion) disk. Omitting the constant term in equation (17) gives 

2E = u> o 2 r 2 + v r 2 + v 0 2 (18) 

Consider now a distribution function of the form 

f(E,J) = C 1 (C 2 - 2E + 2o>J)" 1//2 (19) 

where o> is the constant angular velocity of the disk, constrained so that 0 ^ u> = w Q . 
This requires that 

C 2 - 2E + 2coJ = 0 
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Thus, for r = R, v r = 0, and v 0 = Ra>, 


c 2 - co q 2 R 2 - o) 2 R 2 + 2co 2 R 2 = 0 


or 


c 2 = (co 0 2 - w 2)r2 

Substituting for E, J, and C 2 allows equation (19) to be written as 

- 1/2 

f(E,J) = 0/(^2 - o;2)r 2 - 2E + 2ajjj 


= C 





From equation (15), 


fi( r) = Jj f(E ,J)dv r dv 0 = 27 rRC 1 Ja > 0 2 - w 2 ^ 


ri 

R 2 


( 20 ) 


Therefore, since p.(r) = 11 (0)^1 - the following equation is obtained: 


_ M( 0) 

2ttR^ 0 2 - a; 2 


Ci = 


and 


f(E,J) = f(r,v r ,v 0 ) = 


M(0) 


2ttR 




^ 0 2 - o; 2 )(R 2 - r 2 ) - v r 2 - (v Q - r o>j 


SI -1 / 2 


( 21 ) 


The distribution function (21) is a stationary solution of the Boltzmann equation (12). In 
equation (21) o> = u> 0 corresponds to the cold, violently unstable disk which was previ- 
ously investigated (refs. 2 and 3). For <x> = 0 the disk is nonrotating and purely pres- 
sure supported. 

Toomre (ref. 6) estimated that for a disk with a Gaussian velocity distribution, a 
velocity dispersion equal to or greater than 


a . = 3.36 — 
r,min K 
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should be locally stabilizing for exponentially growing axisymmetric modes. Toomre's 
evaluation of o r min is not easily extended to the present non-Gaussian velocity distri- 
bution given by equation (21). Also, overstabilities may be present since the density of 
stars in phase space is not a decreasing function of epicyclic amplitude (ref. 12). Never- 
theless, equation (18) should represent the local criterion for squelching axisymmetric 
exponentially growing modes (ref. 16). Therefore, 


°r,min 

is used as a parameter in investigating the stability of the disks. The value Q = 1 is 
expected to be a lower bound on the rms velocities needed for stability. Using f^r ,v r ,v^ 
as given by equation ( 21 ) results in ^since v r and \g - ru appear in the same form in 
equation ( 21 ), 0^2 = a r 2 j 


a 0 2 + cr r 2 = 2 cr r 2 


1 

fi(r) 



+ ( v 0 " 


rw ) 2 ] f ( r > v r’ v e) dv r dv e 


= | K 2 • w2 )( r2 - r2 ) 


or 


°r = 


(a) Q 2 - o>2J(r 2 - r2) 


Cl 1/2 


For the uniformly rotating disk, a r m j n is given by equation (5), so that 


Q = 


= 1.69 


Z 1 " ( w / w o)‘ 


( 22 ) 


(23) 


The variation of Q as a function of u> for the uniformly rotating disk is shown in fig- 
ure 24. The value of Q varies from zero for u> = a> 0 to 1.69 for u> = 0; Q is equal 
to 1 for w ~ 0 . 8 co 0 . 
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Results 


The results presented in this section were obtained with a 128 x 128 active mesh 
for the potential calculations. The disks consisted of 100 000 stars and there were 
200 time steps per rotation. 

The case of a cold disk o> = oj Q was previously investigated (ref. 3) and as expected 
the disk was found to be violently unstable. The present investigation considers disks 
with oj = 0.8u> o , a) = 0.6u> o , oj = 0.4<u o , and co = 0 corresponding to initial values of 
Q given by Q = 1.01, Q = 1.35, Q = 1.55, and Q = 1.69, respectively, where Q is 
given by equation (23). It should be emphasized that the model galaxies represented by 
equation (21) do not represent a mass or velocity distribution that one would expect to 
find in nature. For example, for the maximum velocity at a given radius 


v max “ |/K 2 - " 2 )( r2 - r2 ) 

f^r,v r ,Vg|j as given by equation (21) is actually singular, and f^r,v r ,v^ increases with 
increasing v r or vg. However, presently there are no really "good" stationary solu- 
tions of the collisionless Boltzmann equation available for disk galaxies. Some interesting 
solutions for self-gravitating disklike stellar systems are discussed by Miyamoto (ref. 17). 

The evolution of four disks of stars corresponding to equation (21) with 
(a) oj = 0.8w q , (b) oj = 0.6w o , (c) oj = 0Aoj q , and (d) oj = 0 is presented in figure 25. 
Each of the 100 000 stars in the simulation represents 0.84 x 10® solar masses, so that 
the total mass of the disk galaxy is 0.84 x 1C)11 m o (solar masses). The rectangular bor- 
der enclosing the disks represents the active 128 x 128 array of cells used in the cal- 
culations. The initial radius of the disks is 16 kpc. Since the disks become progres- 
sively more stable as the initial velocity dispersion is increased (or oj is decreased), 
the evolution of the more stable systems is investigated for longer times. The times 
shown in figure 25 and subsequent figures are in units of the rotational period of the 
cold (zero-velocity-dispersion) disk T 0 = 2i t/oj q . Figure 25(a) shows that for Q = 1 
(or oj = 0.8aj o j, the system is unstable and within two rotations it has formed a bar- 
shaped structure. After three rotations this structure remains essentially unchanged. 

It should be noted that all small-scale instabilities which occurred in the cold disk (ref. 3) 
have been stabilized. Only the large-scale "bar making" instability is present. A simi- 
lar result is shown in figure 25(b) for Q = 1.35, or oj - 0.6o> . However, the bar struc- 
ture is now much less pronounced. For Q = 1.55, the system is essentially stable. Some 
of the stars near the edge of the disk tend to escape to larger radii. This is to be expected 
since the distribution function f^r,v r ,v^ is singular at the edge and star orbits tend to 
be unstable there. Similar results are obtained for the nonrotating disk shown in fig- 
ure 25(d). Figure 25 indicates that the disk becomes stable for values of Q somewhere 
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between 1.35 and 1.55, or for values of a> between 0.6w o and 0.4a> Q . These results 
are in agreement with a normal-mode analysis of a family of uniformly rotating disks 
which was performed by Agris J. Kalnajs at Tel-Aviv University, Tel-Aviv, Israel. In 
his analysis (as yet unpublished), Kalnajs finds that the simple mode corresponding to 
the bar disturbance becomes unstable for values of u> greater than about 0.508. 

The evolution of the distribution of the radial velocities of the stars as a function of 
star radial distance is shown in figure 26. As can be seen from figures 26(a) and (b), the 
radial velocities of the stars increase rapidly as the system evolves. Also, a large num- 
ber of stars greatly increase their distance from the center of the disk. For c o = O.4o> 0 
and w = 0, the results in figures 26(c) and (d) show that there is little increase in the 
radial velocities, which of course were already large at t = 0. Only a few stars increase 
their radial distance beyond the initial disk radius, especially for the disk in figure 26(d). 
Nevertheless, the distribution of radial velocities at t = 3 in figure 26(d) shows that 
quite subtle changes take place in the velocity distribution near the edge of the disk. 
However, these changes do not appear to affect the structure of the disk appreciably. 

The rectangular border enclosing the velocity distributions in figure 26 extends from 
-350 km/sec to 350 km/sec and from 0 to 30 kpc. 

In order to obtain more quantitative information than can be obtained from figures 25 
and 26, the disk is divided into a number of concentric rings, each 1/2 kpc in width. The 
radial dependence of various parameters averaged azimuthally over each ring is then 
obtained. 

An indication of how hot a disk of stars becomes can be obtained from the evolution 

°T> 

of Q = - — . Figure 27 shows the evolution of the azimuthally averaged Q for the 

°r,min 

four disks. The results in figure 27(a) indicate that the disk becomes rather hot, with 
values of Q in the outer parts of the disk near 5. For the disk with an initial angular 
velocity of w = O.6u> 0 , the value of Q increases from 1.35 to about 2. Smaller increases 
occur for the disk in figure 27(c). Finally, the value of Q in figure 27(d) for the non- 
rotating disk remains nearly constant at 1.69. 

The evolution of the azimuthally averaged density for the four disks is presented in 
figure 28. For the most unstable disk, shown in figure 28(a), the final central density 
increases to a high value given by an approximately exponential density variation. Simi- 
lar results were previously obtained for other violently unstable disk galaxies (ref. 3). 

In figure 28(b) the central density oscillates between 1.4 x 10^ and 2.1 x 10^ M 0 /kpc2, 
reaching a final value of about 1.7 x 108 M Q/ /kpc2 after 4.5 rotations. The changes in the 
density for the disk shown in figure 28(c) were the least pronounced of the four disks. 

For the nonrotating disk in figure 28(d) the central density appears to be oscillating near 
the value 1.4 x 10® M Q /kpc2. When the kinetic energy of the disks is plotted as a function 
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of time, it is found that the kinetic energy initially oscillates nearly sinusoidally about 
the equilibrium value. The initial amplitude of the oscillations is about 5 percent of the 
equilibrium value and the period of the oscillations in all four cases is near (slightly less 
than) the rotational period of the cold balanced disk. For the disks in figures 25(a) and 
(b), the oscillations are strongly damped after the first two oscillations. For the disk in 
figure 25(c), the oscillations are more slowly damped, whereas for the disk in figure 25(d), 
the oscillations in the kinetic energy show no damping during the nine rotations 
investigated. 


CONCLUDING REMARKS 

The results on the evolution of initially balanced disks show that a velocity disper- 
sion given by Toomre's criterion will stabilize a disk against axisymmetric disturbances. 
Also, all fast-growing small-scale disturbances are stabilized. However, all such "sta- 
bilized" disks investigated were found to be unstable against more slowly growing large- 
scale nonaxisymmetric disturbances, and the system consequently assumed a two-arm 
spiral structure. After about three rotations, the spiral structure disappears and the 
central portion of the disks assumes an oval or bar-shaped structure surrounded by a 
hot axisymmetric distribution of stars. In addition to the results presented here, the 
evolution of disks for which the logarithm of the initial density decreased linearly and 
quadratically with radius (i.e., exponential or Gaussian laws) is also investigated. The 
results obtained were similar to those presented here. 

It was possible to generate an axisymmetric stable disk only with considerable dif- 
ficulty. The initial condition for the axisymmetric stable disk was obtained by symme- 
trizing out the bar or oval structure in the central portion of the disk. The resulting 
axisymmetric disk was stable, and the radial variations of all parameters remained as 
they were for the final state of the disk with the central bar structure. 

Attempts at slow cooling of stable axisymmetric disks indicated that whenever 

cr r . 

Q = = — (the ratio of the velocity dispersion to the minimum velocity dispersion given 

a r,min 

by Toomre's stability criterion) reached about 2, further cooling would only cause the 
occurrence of collective instabilities which would heat up the disk as fast as it was being 
cooled. 

An interesting end result for all the disks of stars (which evolved from an initially 
unstable state) investigated so far is that the final distribution in the radial direction for 
the "disk population of stars" is closely approximated by an exponential variation of den- 
sity. This result may be significant since it agrees with observational evidence which 
indicates that the luminosity in the outer regions of many spiral and SO galaxies seems 
also to decrease exponentially with radius. 
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The attempt to induce spiral structure by means of the passage of a companion 
galaxy was only partly successful. Because of the rather large velocity dispersion of the 
disturbed galaxy the induced spiral structure was weak. 

The investigation of a disk of stars which is a stationary solution of the collision- 
less Boltzmann equation showed that the disk became unstable at a value of the velocity 
dispersion which agreed with the predictions made by a normal- mode analysis. For 
rotational velocities a> of 0.8w o and 0.6co o (where oj q is the rotational velocity of the 
cold, zero-velocity-dispersion disk) corresponding to Q = 1 and 1.35, respectively, the 
disk formed a barlike structure. For u> = 0.4cu o and u> = 0 (or Q = 1.55 and 
Q = 1.69), the disk was stable against the bar-forming mode; however, for these two 
cases the disks sustained what appear to be natural oscillations, or pulsations, with a 
period near the rotational period of the cold balanced disk. The pulsations were espe- 
cially pronounced for the nonrotating disk, where they did not show a decrease in ampli- 
tude for the nine rotations investigated. Since the initial conditions were generated by 
means of a pseudorandom-number generator, small oscillations are to be expected. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 10, 1972. 
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APPENDIX 


COMPUTER MODEL 
Discretization Parameters 

The model for the disk galaxy consists of a large number of representative stars 
(here 100 000) that are confined to move in the galactic disk. An N x N (here 
256 X 256) array of cells is superposed over the plane of the disk for the purpose of cal- 
culating the gravitational potential. At the center of each cell a mass density is defined 
which is given by the number of stars in that cell. The mass -density distribution is used 
to obtain the gravitational field at the center of each cell. From the gravitational field 
the force acting at the position of a star is calculated by means of a bilinear interpola- 
tion among the fields of the four cell centers surrounding the star. Newton's equations 
of motion are then used to advance the position and velocity of each star by a small time 
step. Typically, there are 200 time steps per "galactic rotation." If a star should leave 
the N x N array of cells, it is still included in the calculations by approximating the 
force acting on the star. With the exception of the improved potential solver described 
in this appendix, the model used is the same as that previously described (ref. 3). 

The effects of varying the number of stars and various other discretization param- 
eters were investigated by Hohl and Hockney (ref. 2). The effect of binary collisions for 
the model has been estimated (ref. 3) to be such that the collision time equals about 
100 "galactic rotations." The time history of the velocity dispersion of stars near var- 
ious radii presented in figure 14 also indicates that collisional or thermalization effects 
are unimportant for the time period used for the present calculations. 

Potential Calculation 

The scaled gravitational potential at the center of cell (x,y) is defined by the double 
summation over the two-dimensional array of cells 


where 


N-l N-l 

^x,y = Z 2 

i=0 j=0 


Hi, j =(i 2 + j 2 )“ 1/2 


H 0,0 = 1 


(Al) 


(i + j * 0) 
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APPENDIX — Continued 


and py is the mass density in cell (i,j). The double summation is evaluated by the 
convolution method using fast Fourier transforms (ref. 2). That is, the Fourier trans- 
form of the potential equals the product of the Fourier transforms of p and H 




(A2) 


The gravitational potential <p v is obtained by taking the inverse Fourier transform of 

equation (A2). Rather than using a complex Fourier series, a real expansion was chosen. 

For example, the Fourier transform of the density p is given by 

x,y ** 


5 M= Y Y 

y=0 x=0 


(O^k, is n ) 


N-1N-1 

I I c(x)p x>y cos(^)sin[7r(Z-n)|] 
y=0 x=0 


(0 S k g n) 
(n < l < N) 


N-l N-l 

I I c(y)p X5y sin[7r(k-n)|]cos(^.) 
y=0 x=0 


> (A3) 


(n < k < N) 
(O^Sn) 


N-l N-l 

^k,Z = ^ Z 4 x,y sin [^ k - ^llsinjua - n)|l (n < k, l < N) 

y=0 x=0 


where 


c(x) = 1 j\[2 if x 


0 or x = n 


c(x) = 1, otherwise 


The symbol n defines the n x n active array and N = 2n defines the larger array 
over which the Fourier transform must be taken so that the potential for an isolated 
disk galaxy is obtained. Note that the density may be nonzero only in the smaller 
n x n array. Because of the symmetry of H x>y , the Fourier transform H k l can 
be obtained by a finite cosine transform: 
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APPENDIX — Continued 


n n 


fi k ,1 = 22 c2 ( x ) c2 (y) H x,y cos (^) cos (?) 

y=0 x=0 


(0 g k, l Z n) (A4) 


and 


H k+n,Z " H k,Z+n H k+n,Z+n = H k,Z 


The next step in obtaining the potential is to multiply ^ by ^ to obtain 
K,l = ^k,Z**k,Z 

The gravitational potential for an isolated galaxy correctly defined over the n x n array 
is obtained by the Fourier synthesis 




N-l 

cos (f kx ) + E 

k=n+l 


\ i sin [l (k - 



+ 


N-l n 

II 


i=n+l k=0 


\i cos {^) + 


N-l 

I 


k=n+l 





(A6) 


A FORTRAN listing of the computer program actually used to obtain the potential by use 
of an N x n array of cells is given on the following page. The variable I2A defines 
the size of the rectangular array used for the potential calculations. When the subroutine 
GETPHI is called, RHO(I,J) contains the mass density and GETPHI places the values of 
the corresponding gravitational potential in RHO(I,J). The subroutine FTRANS(I,I2B) was 
written by Hockney (ref. 18) and it performs a finite Fourier analysis or synthesis on the 
COMMON input array Z and places the result in the COMMON output array Y. The 
subroutine performs a cosine analysis for I = 2, a periodic analysis for 1 = 3, and a 
periodic synthesis for 1 = 4. The subroutine GETSET(I,I2B) initializes FTRANS and 
is called every time the arguments of FTRANS(I,I2B) are changed. The Fourier trans- 
form ^ is calculated on an (n+1) x (n+1) array only the first time that the subroutine 
is called and is kept in storage for subsequent use. 
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APPENDIX - Continued 


SUBROUTINE FOR CALCULATING THE GRAVITATIONAL POTENTIAL 


SUBROUTINE GETPHI 

COMMON Z(257) ,Y(257) ,RH0(256,128) , I2A, ITEST 
DIMENSION H(1 29 ,T 29) 

IF(ITEST.EQ.O) GO TO 10 

ITEST=0 

I2B= I2A- 1 

N=2**I2A 

N02=N/2 

N21=N02+1 

RNI=1./(N*N) 

DO 1 J=1 ,N21 
DO 1 I=1',N21 

I F ( I . EQ. 1 . AND. J. EQ.l ) GO TO 1 

H(I ,J)=RNI/SQRT( (1-1 . )*(I-1 . )+(J-l . )*(J-1 . )) 

1 CONTINUE 
H(1 >1 )=RNI 

CALL GETSET(2,I2B) 

DO 2 J=I ,N21 
DO 3 1=1 ,N21 

3 Z( I )=H ( I ,J) 

CALL FTRANS(2 , I2B) 

DO 4 1=1 ,N21 

4 H { I ,J )=Y ( I ) 

2 CONTINUE 

DO 5 1=1 ,N21 
DO 6 J=1 ,N21 

6 Z(J )=H( I ,J) 

CALL FTRANS(2,I2B) 

DO 7 J=1 ,N21 

7 H ( I ,J)=Y(J) 

5 CONTINUE 
10 CONTINUE 

CALL GETSET(3,I2A) 

DO 11 J=1 ,N02 
DO 8 1=1 ,N 

8 Z ( I )=RHO(I ,J) 

CALL FTRANS ( 3 , I 2A ) 


DO 9 1=1 ,N 
9 RHO(I ,U )=Y ( I ) 

11 CONTINUE 

DO 12 1=1 ,N 
DO 13 J=1 ,N02 
Z(J)=RHO(I,J) 

13 Z(J+N02)=0. 

CALL GETSET(3,I2A) 

CALL FTRANS(3,I2A) 

I F ( I . GT . N21 ) GO TO 14 
DO 15 J=2 ,N02 
Z(J )=Y (U )*H(I ,0 ) 

15 Z(J+N02)«Y(J+N02)*H(I,U) 

Z(1 )=Y(1 )*H(I ,1 ) 

Z(N21 )=Y(N21 )*H(I ,N21 ) 

GO TO 16 

14 DO 17 J=2 ,N02 
Z(J)=Y(J)*H(I-N02,J) 

17 Z(J+N02)=Y (J+N02)*H(I-N02 ,J ) 
Z ( 1 )=Y(1 )*H( I-N02 ,1 ) 

Z(N21 )=Y(N21 )*H(I-N02,N21 ) 

16 CONTINUE 

CALL GETSET(4,I2A) 

CALL FTRANS (4, I2A) 

DO 18 J=1 ,N02 

18 RHO(I,J)=Y(J) 

12 CONTINUE 

DO 19 J=1 ,N02 
DO 20 1=1 ,N 

20 Z(I)=RHO(I,J) 

CALL FTRANS (4, I 2A) 

DO 21 1 = 1 ,N21 

21 RHO( I ,J )=Y( I) 

19 CONTINUE 
RETURN 
END 


Next the Fourier transform of in the x-direction is obtained on the 

N x n array, that is, for 0 ^ x ^ N - 1 and 0 in, Since fx x y is nonzero 
only over the n x n array, the components of the Fourier transform of ^x,y *- n 
the x-direction will be zero for n < y < N. Therefore, by use of the one-dimensional 
arrays Y and Z, one can perform the Fourier transform in the y-direction, multiply 
the result by and take the inverse Fourier transform in the y-direction. The 

result is placed in the N x n RHO(I,J) array for 0 = y ^ n - 1 and 0 ^ x ^ n - 1 
with the values for n < y < N discarded. The final step is to perform the inverse 
Fourier transform in the x-direction for 0 In - 1. This procedure gives the cor- 
rect gravitational potential for an isolated disk galaxy over the n x n array. 

Table I gives the measured time for calculating the gravitational potential with the 
program listed above. Also shown are the times required for the previous potential sol- 
ver which required a larger N x N array. A listing of that program is given in 
reference 3. 
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APPENDIX - Concluded 


TABLE I.- COMPUTER TIME REQUIRED TO OBTAIN 
THE GRAVITATIONAL POTENTIAL 


CDC 6600 CPU time, sec, for — 


Active 
n x n 
mesh 

Present potential solver 
using 2(n x n) storage 

Previous potential solver 
using 4(n x n) storage 


H calculated 

H given 

H calculated 

H given 

16 X 16 

0.164 

0.126 

0.182 

0.144 

32 x 32 

.618 

.476 

.660 

.526 

64 x 64 

2.440 

1.892 

2.626 

2.086 

128 x 128 

10.000 

7.740 

10.794 

8.530 
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Computer simulation 


M 81, Sb-type spiral galaxy 


L-72-125 

Figure 1.- Comparison of four computer- generated galaxies with photographs of actual galaxies. 






Computer simulation 


NGC-175, SBab-type barred spiral galaxy 

L- 72-126 


Figure 1.- Continued. 



Computer simulation 


M 101, Sc-type regular spiral galaxy 

L- 72-127 


Figure 1.- Continued. 
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Computer simulation 



NGC-1073, SBc-type barred spiral galaxy 

L-72-128 


Figure 1.- Concluded. 





t = 0 t = 1.6 t = 3.2 



t = 4.8 t = 6.4 t = 8.0 


Figure 2.- Axisymmetric evolution of an initially balanced uniformly rotating 
disk of 100 000 stars. The stars have an initial velocity dispersion given 
by Toomre's criterion and move under a purely radial gravitational field. 
Time in this and all subsequent figures is given in units of the rotational 
period of the cold balanced disk. 
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t = 6.0 


t = 8.0 


t = 10.0 


Figure 5.- Unconstrained evolution of the initially balanced uniformly 
rotating disk of 100 000 stars. The stars have an initial velocity 
dispersion given by Toomre's criterion. 
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Figure 6.- Evolution of the azimuthal variation of the star density at three 
radii r (in kpc) for the disk shown in figure 5. 
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Figure 7.- Variation of the radial velocity dispersion as a function of 


radius for the disk shown in figure 5. 
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Figure 8.- Dependence of Q = ^r/ a r,min on rac *iiis f° r the disk shown in figure 5. 
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Figure 9.- Evolution of the mass density plotted as a function of radius. The 
density is given in units of solar masses per kpc2. (Note that each one of 
the 100 000 simulation stars has a mass of 0.84 x 10®M G .) 
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• Density across bar 
■ Density along bar 
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Figure 11.- Density variation along and across the bar at t = 8.0. 
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Figure 14.- Evolution in time of the rms velocity at four radii for the 
stable axisymmetric disk shown in figure 12. 



r, kpc 

Figure 15.- Comparison of the radial variations of the angular velocity co = ^K r /r, 
the azimuthal velocity ru, and the mean azimuthal velocity of the stars (V q) . 
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Figure 27.- Evolution of Q = u r /<V m in for the four disks shown in figure 25 for 
(a) w = 0.8w o ; (b) w = 0.6co 0 ; (c) o> = 0.4w o ; and (d) w = 0. 
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Figure 28.- Evolution of the azimuthally averaged density for the four disks shown in figure 25 for 
(a) u) = O.8o> 0 ; (b) co = O.6u> 0 ; (c) co = 0.4co o ; and (d) co = 0. 



NATIONAL AERONAUTICS AND SPACE AD M I ST RAT 1 0 N 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 
PENALTY FOR PRIVATE USE $300 


FIRST CLASS MAIL 


postage and fees paid 

NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 



025 001 Cl U 30 720317 SOOQD^nc 
DEPT OF THE air FORCE S00903DS 


AF WEAPONS LAB (AFSC) 

TECH LIBRARY/WLOL/ 

ATTN: E LOU BOWMAN, CHIEF 
KIKTL AND A FB NM 87117 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual ) Do Not Return 


'The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atyno sphere and space. The Administration 
shall provide for the widest practicable and appropriate disse??iination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



; & V 




I 




